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Abstract. We extend the curvaton scenario presented by Erickcek et al. [1,2], to explain how the 
even-odd multipole asymmetry of the Cosmic Microwave Background (CMB) (also called parity 
asymmetry, [3, 4]) and power anisotropics can be generated by the curvaton field, which acts as an 
extra component to the spectrum of adiabatic perturbations in the inflationary epoch. Our work pro- 
vides a possible cosmic explanation to the CMB large-scale asymmetry problems besides systematics 
and unknown residuals. 
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1 Introduction 

The Cosmological Principle states that the Universe should be isotropic and homogeneous on scales 
above 100 Mpc, which is widely accepted as a basic principle of most cosmological scenarios. This 
principle can be experimentally tested by galaxy surveys and CMB observations. The SDSS ex- 
periment result indicates that the galaxy distribution becomes isotropic and homogeneous at large 
scales [6-12], which supports the Cosmological Principle well. 

However, there are still open questions from the CMB observations-so-called anomalies. Among 
them is the non-Gaussian cold spot [17-19], the missing power in the quadrupole in all WMAP 
data [5, 13-16] (see however [5]), even-odd multipole power asymmetry (also called parity asym- 
metry [3, 4]), alignment of multipole components and axis-of-evil [20-23], north-south asymme- 
try [24-28] and so on. Certainly, there are attempts in the listed literatures and others to attribute 
those anomalies to systematics or unrealized Galactic/Ecliptic emission, or simply the cosmic vari- 
ance. The systematics can be reduced or eliminated by cross check between different and independent 
CMB observation experiments, like WMAP and Planck, thus we should wait patiently for the Planck 
data release for a conclusion regarding the systematics issues. As for the cosmic variance, it could 
explain many anomalies but only with very low probability. There is an even more fundamental ques- 
tion: whether these anomalies have a common origin or are statistically independent. If we believe 
that all the anomalies originate from the same source, it would be even more important to discover 
its origin. In the opposite case, where all the anomalies are statistically independent, the problem 
is how one peculiar realization of the random field can contain all these anomalies. Although the 
anomalies have been observed in the temperature data, they could also prove to be main sources of 
contamination in polarization data. Therefore, the understanding of their origin is potentially crucial 
for investigating the E-mode and B-mode or even the primordial non-Gaussianities and gravitational 
waves from inflation. 

Therefore, if we believe that the anomalies of the CMB discussed in the literature indicate some 
constraint on the Cosmological Principle, we should find an explanation through cosmological theory. 
In the standard inflationary scenario, the large-scale structure is generated by the initial perturbations 
due to quantum fluctuations of the inflation field. However, if we further consider the possibility 
that the standard inflation field is not the only field in the inflation stage (the inflation field is still 
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dominating), then by adding some extra components to the nearly scale-invariant spectrum we can 
introduce a seed of asymmetry to the theoretical expectation, not just to a specific realization by a 
particular observer, which is the main purpose of this work. 

In the curvaton scenario [29, 30], we can see that additional non-isotropic perturbations can be 
generated by the curvaton field, and consequently cause even-odd multipole power asymmetry in the 
power spectrum. In the curvaton scenario, the curvaton field (a) is supposed to have negligible energy 
density compared to the inflation field. It is also non-interacting with the inflation field, and thus its 
initial value a* is kept during inflation, and its quantum fluctuation {Sa)rms = ^fm//(27r) (Hinf is 
the Hubble parameter during inflation) contributes part or all of the primordial perturbations [31-34]. 
If the curvaton potential is V{a) = (l/2)m^cr^ with rricr <C i^m/ (^a is the mass of the curvaton), 
then after inflation (where ma — H) the curvaton will oscillate and decay into radiation and will 
interact with matter. The sequence of curvaton decaying and decoupling of particle species gives 
different curvaton interacting scenarios, like curvaton- dark matter interacting [2]. 

The curvaton scenario discussed in this work is an extension of [1, 2]. Whereas Erickcek et al. 
focus on super horizon perturbations (the wave length of the perturbation they consider is very large), 
we have discovered that if the wave length of the curvaton perturbation is comparable to or smaller 
than the horizon, then the model can be used for explanation of some of the CMB anomalies. 

The outline of this paper is the following. We present the extended model in section 2. In section 
3, we apply this model to the WMAP data to see if it can, at least partly, explain the power spectrum 
parity asymmetry and the temperature space anomalies. In section 4 we show how a plane wave 
component can affect the CMB power spectrum. In the end, a brief discussion is given in section 5. 



2 Extended model based on the curvaton scenario 



Based on the curvaton scenario, we have constructed a model with only three parameters to see if it 
can generate some of the observed anomalous features of the CMB, in particular power asymmetry 
in the power spectrum. The model is presented below. 

Following to [4], if there is a primordial perturbation in Fourier space ^^(k), then the low mul- 
tipole (2 < / < 30) spherical harmonic decomposition coefficients {aim) are connected to $(k) 
through 

^*(k)TKfc)y4(k), (2.1) 

where Ti{k) is the radiation transfer function. For odd multipoles, / = 2n + 1 (n = 0, 1, 2, ..), 

aim = -^^J^ / dkk" j d9^sm9^ x J #krKA:)l^;,(k)Im[*(k)], (2.2) 



and for the even multipoles, / = 2n, 
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J dkk"^ J d9u sin Oux J d<j)i,Ti{k)Yi*^{k)Re[^{k)]. (2.3) 



From Eq.(2.2) and (2.3), we can see that e.g. odd-parity preference might be produced, provided that 

|Re[*(k)]| « |Im[*(k)]| {k<22M, (2.4) 
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where tjq is the present conformal time. As is seen from 2.4 the phases of metric perturbations 



— arctan 



Im(^(k)) 
Re(^(k)) 



) have to be locaHzed in the vicinity of ^ oc 7r/2; 37r/2, at least for the range 

k < 20/?7o to 30/770. This shows the possibiUty of generating even-odd parity asymmetry from 
specific primordial perturbations. 

The squeezed space of phases indicates that for spatial scales x > 4 Gpc [4] the homogeneity 
and isotropy of the perturbations is abnormal. Namely, parity arguments of the CMB leads to the 
parity asymmetry of the metric perturbations, ^^(x) ^ — ^^(— x). Let us assume that the origin of 
those anomalies can be associated with unusual properties of the curvaton field [1-3]. The potential 
perturbation at decoupling due to a curvaton field perturbation is given in [1], using a real space form 
"^{rdec^x), as: 
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(2.5) 



where a is the homogeneous curvaton background, and a{x) = a + Sa(x). R = Pa/ Ptotai is the 
fraction of curvaton in the total energy density just before curvaton decay. The curvaton decay is 
assumed to be early enough so that 2^ —{2R/9)Spcr/ Pa- 

We set the time-dependent coefficient of [2(^) + (^)^] as ^(r) and suppose ^ = r sin(fc • 
X -\- S) (sinusoidal fluctuation for curvaton perturbation, see also Sec. 4 of [2], and S is the phase), 
then we have the spatial distribution of the real space potential as: 



(2.6) 



^(r,f) = '0(T)[2r sin(fc • f + 5) + r^sin^(fc • x + S)]. 

The low order CMB power spectrum consist of Sachs-Wolfe (SW) and Integrated Sachs-Wolfe 
(ISW) effects. According to [1] the induced SW effect is = ^(^(iec)/3. Therefore, we 

only have to calculate the ISW effect. According to Equation 16 of [1], the ISW effect is given by: 

= 2/' ■"l^{».g.-'lx.-xWl»}^_ ,,,, 

ISW Jadec 

where x{o) = Ho[T{a) — Tdec] and xo = = 1) = H^Xdec If we assume that r is constant this 
gives: 
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(2.8) 



Combining the SW and ISW effect, and letting k • Xdec = ^7r[l — cos(0)] = qou (q is the wave 
number, 9 is the polar coordinate, and we choose the system of coordinates so that k is oriented along 
-Z), we have gotten the CMB fluctuations due to curvaton perturbations as: 



AT ■ 



= 4r\l>c[sin(g<^ + 5) + - sin^(g<^ + 5)], 



(2.9) 



where ^'c = ^{a) 1^^^^ +^(T(iec)/3 is a constant and is related only to the overall amplitude. 

Now it is clear that in our single-wave model, the structure of CMB fluctuations due to curvaton 
perturbations are determined by only three parameters: q, r and S. At the current stage, we consider 
only the structural term of Equation 2.9: 



AT 



(n) 



oc sm{q(jj + 5) + - sm'^{q(jj + 5), 



(2.10) 



and let the amplitude from equation 2.9 be a free parameter. From this equation we can get the 
curvaton component power spectrum. 



-3- 



3 Implementation of the curvaton model 



The curvaton model can now be implemented and compared to CMB data. The model is determined 
by three parameters: the wave number q, the curvaton fluctuation strength r and the initial phase 
5. Note that changing 5 is very similar to choosing a special spatial position of a particular observer. 
Firstly, in section 3. 1 we determine the model parameters by fitting the WMAP power spectrum. Then 
in section 3.2 we proceed to find the most optimal orientation of the model based on the WMAP data. 

3.1 Determining model parameters by fitting the WMAP power spectrum 

We apply the model to the WMAP CMB power spectrum. The best fit ACDM power spectrum does 
not have power asymmetry, but the observed WMAP CMB power spectrum does. Thus, assuming 
that the WMAP CMB power spectrum is a combination of ACDM and an extra component due 
to the curvaton field, the power spectrum of this extra component {Cextra) should display power 
asymmetry. Therefore we fit our model to the difference between the ACDM best fit power spectrum 
and the observed WMAP CMB power spectrum. 

With each parameter set (g, r, 5), we calculate the CMB temperature distribution according to 
Equation 2.10 as well as the CMB power spectrum, and then linearly fit the derived power spectrum 
to Cextra to determine the constant (Equation 2.9). The statistic of fitting is recorded for 
this parameter set, and the best guess of (g,r, 5) is determined by the minimal x^- The resulting 
CMB power spectrum is given in the top left of Fig. 1. It seems as though the characteristic power 
asymmetry structure has been faithfully produced by the model, and only the / = 2 (quadrupole) 
component is not particularly well fitted. To quantitatively estimate whether parity asymmetry is 
actually produced, we use an estimator g{l) (see [3, 4]). The estimator is defined as: 

Y!iZTi{i + i)c^' 

where C+(/) = C{1) cos^(y) (the power spectrum for all even /) and C~{1) = C{1) sm^{'^) (the 
power spectrum for all odd /). Note that we take the sum from / = 3 because of the poor fit of the 
quadrupole. In Fig. 1 (bottom) we see that the parity asymmetry {g{l) < 1) is indeed reproduced for 
the model for low multipoles, as expected. 

We also give a contour plot of the likelihood of fitting in the parameter space. The likelihood 
of fitting can be calculated as Cfu = P{x^ > x}it)' have tested that the model CMB power 
spectrum is not sensitive to r and thus we choose to fix r at its best guess value, r = 2.6, to plot a 
2D-contour of ln{Cfit) over q and S. This is given in the right panel of Fig. 1. We can see that there 
is a double-peak structure along the S axis, separated at about O.Ttt. We have confirmed that these 
two peaks give very similar resulting power spectra. It is not strange to see a double-peak separated 
by O.Ttt, because all large-scale perturbations are more or less spatially periodic. 

Although the fitting in Fig. 1 looks nice, we must be careful about concluding that the entire 
large-scale CMB asymmetry is generated by our model. For now, we can only say that part of the 
large scale power asymmetry can be explained by our toy-model. 

3.2 Determining the most optimal orientation of the model 

Since the temperature fluctuations caused by plane wave curvaton perturbations are rotationally sym- 
metric around the wave vector k (see Fig. 4 for an example), its spherical harmonic components ai^ 
should not have the same strength at different m. Especially, if k has the same direction to ±Z-axis, 
then aim — for all m 7^ components. Therefore, if the model in this work is real, then the 



-4- 




3 4 5 6 7 8 9 10 11 12 
L 



Figure 1. Application of our model to WMAP data: Top left: The best fit model power spectrum (dotted) 
and the difference between the ACDM best fit power spectrum and the observed WMAP 7-yr power spectrum 
(solid). Top right: A 3-level contour plot of ln{Cfit) at 0.5, 0.75, 0.96 (black, blue, red), where ln{Cfit) is 
normalized to (0, 1). Bottom: Parity parameter for observed WMAP 7-yr power spectrum (solid line) and for 
the curvaton model (dashed line). 



orientation that minimizes the m 7^ spherical components for the real CMB data is very likely the 
orientation of the model. According to Fig. 1, we see that for / = 3 ^ 7 components, the parity asym- 
metry is most significant, thus if we work in this range, it may increase the accuracy of determining 
the orientation of the model. 

Our approach is like this: First we rotate the WMAP 7-year ILC map around the Y-axis (Galac- 
tic plane) to find an angle 9i that minimizes X^m/o l^^mP for each / in range / = 3 ^ 7 respectively. 
The average value {9i) tells us the latitude of the orientation in the Galactic coordinate system, which 
is —53°. Interestingly, as discovered by [36] the preferred axis of the WMAP quadrupole (/ = 2) 
and octupole (/ = 3) both point to (/, b) ^ (110°, 60°) in Virgo. Since the preferred axis does not 
distinguish between h and —h (b = ±60°), we see that the axis we have found is only 7° away in 
latitude from the well known "axis-of-evil". Moreover, as discovered by [17, 18], the well known 
non-Gaussian cold spot at (/, b) = (209, —57) is also only 4° away in latitude from our axis here. 
Thus our model may play an important role not only in the parity asymmetry problem, but also in 
other well known large-scale CMB anomalies, and perhaps even be connected to some CMB non- 
Gaussianities. 

The standard deviation of 9i in range / = 3 ^ 7 is ctq = 15.4°. Such a small value means 
that these harmonic components have clustered orientations. To confirm that ctq is really small, we 
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Figure 2. Overlapping of belts corresponding to \0 — 7t/2\ = (Oi) for I = 3 ~ 7, plotted in the Galactic 
coordinate system. The two poles of the axis with the strongest overlap is marked out by "lA, IB", and the 
two poles of the secondary axis is marked out by "2A, 2B". 



did a test with 5000 simulations, and for each i-th simulation, the standard deviation of 9i in range 
/ = 3 ^ 7 was calculated as a^. We only got 41 out of 5000 simulations that had < ctq. This 
means that for the WMAP data the clustering of 9i in the range / = 3 ^ 7 is significant at a level of 
99.2%. 

After determining the latitude of the orientation, we change the coordinate system of the ILC 
map by rotating the Z-axis to all 192 directions defined by the HEALPix resolution Ngide = 4 [35]. 
For each new coordinate system, we do the same as we did above and get {9i) for this coordinate 
system. If the harmonic components are sufficiently clustered (i.e. if ctq < 20°) we draw a belt at 
{9i) with a width of 12° for this coordinate system. The overlapping of these belts are shown in 
Fig. 2 (all turned to Galactic coordinates). The hottest spots give the possible orientations of the 
model. The two poles of the strongest orientations are marked out by "lA, IB". The coordinate 
of the axis is (/, b) = (189°, -55°). We can also see an orientation at (/, b) = (346°, -50°) from 
Fig. 2, whose poles are marked out by "2A, 2B". To determine the most optimal orientation of the 
model we calculated the correlation coefficients (for / = 3 ^ 7) between the multipole components 
of our model and the ILC, when the model was rotated to directions "IB" (axis-1) and "2A" (axis-2) 
respectively. As is seen from Table. 1 we have determined that axis-1 is the most optimal orientation. 

3.2.1 Similarity between harmonic components 

We show the similarity between the large-scales components of the model (rotated to axis-1) and the 
WMAP ILC map in Fig. 3. Let us take the / = 5 component (fourth row) as an example: both model 
and the ILC have one cold spot located at —90° to —60°. We also see a band of cold spots between 
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I 


2 


3 


4 


5 


6 


7 


Axis-1 
Axis-2 


-0.13 
0.15 


0.55 
0.15 


-0.16 
0.09 


0.60 
-0.08 


0.59 
-0.01 


-0.05 
0.14 



Table 1. The model-to-ILC correlation coefficients between their multipole components when the model 
temperature map is rotated to the "IB" (axis-1) and "2A" (axis-2) directions. 

-30° and 0° and between 30° to 60°. There is one hot spot located at 60° to 90°, and hot spots in the 
band between 0° and 30°, and the band between —60° and —30°. We can conclude that our rotated 
model fits the structure of the ILC well. 

3.3 Stability of the model to contaminations 

We have tested the stability of our method by using an input map derived from Equation (2.9) with 
g = 3, r = 0, 5 = 0, A: = (1,0,0) (neglecting the coefficient 4r^^c)- By using the same method 
mentioned in the previous section, the resulting axis is found to be (/, 6) = (4.5°, 6°), well consistent 
with expectation, and validates the stability of our method. Then we added two contaminations to test 
the stability of this method: The contaminations are similar to the source but with 1/10 strength and 
different wave vector directions: k = (0, 0, 1), \/3/3(l, 1, 1) respectively. With the same approach, 
we got a resulting axis at (6, /) = (10°, 357°), also close to expectation. Thus we have shown that 
our approach is insensitive to weak contaminations with "misleading" axes (e.g. due to other sources 
of asymmetry). 

4 The effect on the power spectrum of a plane wave component 

Here we show an example from simulation, in which we can clearly see how the power spectrum 
power asymmetry can be generated from a plane wave component due to the curvaton model dis- 
cussed above. First we generated a simulated CMB map from a ACDM power spectrum. Then we 
generated a map of a plane wave component with the same parameters as shown in Fig. 1 (the direc- 
tion of the plane wave is not important for a simulated map, so we choose the Ecliptic north/south 
poles as the direction). The summation of them resembles the "ACDM + curvaton" scenario. We 
then calculate the power spectrum of the original map, do the summation and calculate the total 
power spectrum, and plot them together with the WMAP CMB low-/ power spectrum in Fig. 4. We 
can see that the odd-even multipole power asymmetry of the WMAP low-/ power spectrum can be 
very well reproduced in this way (with the exception of the quadrupole (/ = 2)). We have thus shown 
that the model presented in this work can reasonably account for the WMAP low-/ power asymmetry. 

Certainly, with more simulations we can also see cases in which the combined power spectrum 
(blue line in Fig. 4) isn't similar to the real data (red line in Fig. 4). The reason is simple: the CMB 
and curvaton components can have different directions and phases, thus the summation of them can 
make the power spectrum either higher or lower. This fact makes the problem much more complex. 
However, with more simulations we can see that with the curvaton component presented in Fig. 1, 
the probability of getting similar result to real data will increase. 

The power spectrum similarity is evaluated by the cross correlation coefficient between the 
power spectrum for WMAP data and simulations in the range / = 4 ^ 12 for 10,000 simulations: 

C4_i2 = Corr(Cf ^, cf^^^), / = 4 - 12. (4.1) 

If the simulations are pure ACDM (no curvaton component is added) only 2.6% ± 0.16 of the simu- 
lations have a C4_i2 > 0.6. For the simulations with ACDM + curvaton component 18.4% ± 0.43 of 
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Figure 3. Maps of / = 2 ~ 7 (top to bottom). (Left:) components of the model with the best-fit parameters 
shown in Fig. 1 rotated to the direction of "IB" (axis-1). (Middle:) the WMAP 7-year ILC map. (Right:) 
the model components rotated to the direction of "2 A" (axis- 2). The lines indicate latitudes in the coordinate 
system of the model (the lines on the ILC maps correspond to the coordinate system of axis-1). 
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Figure 4. Top left: The original simulated CMB map with ACDM power spectrum. Top right: The plane 
wave component due to a curvaton field. Bottom Left: The summation of the two top panels. Bottom right: 
The harmonic power spectrum of: theoretical ACDM (Black dash), simulated map without plane wave (black 
solid), with plane wave (blue), and WMAP low-/ power spectrum (red). 

the simulations have a C4-12 > 0.6. This fact supports the curvaton scenario presented in this work 
quite well. 

5 Discussions 

In this work we introduce a model based on the curvaton scenario, which has only three parameters, 
and try to apply it to the WMAP data. This model is an extension of [1, 2], and the main difference is: 
We have discovered that if the wave length of the curvaton perturbation is comparable to or smaller 
than the horizon (q > 1), then the model can be used to explain parity asymmetry and probably 
more asymmetry problems. Our results show that such a simple model can give a very well fit to the 
difference of the CMB power spectrum between ACDM expectation and experimental detection. The 
spatial structure can also be well fitted, especially for / = 5. This tells us that at least part of the CMB 
large-scale asymmetry can be attributed to an extra component to the inflation field, which provides 
a possible cosmic explanation to the CMB asymmetry problems. However, these morphological 
features could also be mimicked by some combination of foreground residuals, as discussed in [37]. 
In this work our goal has been to propose a theoretical model based on the curvaton scenario, and to 
see which constraints we have to apply to it in order to explain the observed parity asymmetry in the 
CMB. 

A recent paper has placed limits on the semi-classical fluctuations in the primordial Universe [38]. 
Although the fluctuation amplitude of our model is not discussed here, we notice that they have ob- 
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tained nearly same direction to us in their Fig. 3 (compare to the "lA" direction in our Fig. 2). 
Therefore, the limits on fluctuation amplitude discussed in their paper may also apply to our model. 
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